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Aims. A magneto-hydrostatic model is constructed with spectropolarimetric properties close to those of solar photo- 
spheric magnetic bright points. 

Methods. Results of solar radiative magneto-convection simulations are used to produce the spatial structure of the 
vertical component of the magnetic field. The horizontal component of magnetic field is reconstructed using the self- 
similarity condition, while the magneto-hydrostatic equilibrium condition is applied to the standard photospheric model 
with the magnetic field embedded. Partial ionisation processes are found to be necessary for reconstructing the correct 
temperature structure of the model. 

Results. The structures obtained are in good agreement with observational data. By combining the realistic structure 
of the magnetic field with the temperature structure of the quiet solar photosphere, the continuum formation level 
above the equipartition layer can be found. Preliminary results are shown of wave propagation through this magnetic 
structure. The observational consequences of the oscillations are examined in continuum intensity and in the Fe I 6302A 
magnetically sensitive line. 
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1. Introduction 
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The study of magnetic elements at very small scales is one of 
, the most important topics in solar physics. Magnetic Bright 
\ Points (MBPs) are ubiquitous in the solar photosphere. 
They have small diameters, typically less than 300 km, and 
are found in the intergranular lanes. MBPs correspond to 
areas of kilogauss fields, are best observed in G-band disk 
centre images and are numerous in active regions or near 
' sunspots. They are formed by a complex process involving 
the interaction of the magnetic field with the convectively 
unstable hot plasma. The physical proc esses associated 
\ with t h eir formation have b ee n outlined in ISchiissler et al.l 
(|2003l ): IShelvag et all ([2004D : ICarlsson et all ([20041 ). using 
■ forward modelling of radiative magneto-convection in the 
solar photosphere and upper convection zone. 

Despite the overall success of photospheric and sub- 
photospheric radiative magneto-convective models to re- 
produce many of the observational properties of solar radia- 
tion, we still do not fully understand the physical processes 
involved in the strongly magnetised photospheric plasma. 
In particular, it is difficult to use the results of these sim- 
ulations for studies of acoustic wave propagation through 
the solar atmosphere and interior. Strong convective mo- 
tions of the photospheric plasma can hide the signatures 
of acoustic waves, making them a difficult subject in both 
numerical and observational investigations. 
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The development of new methods 
the properties of solar plasma using 
have been followed by the successful 
the magneto-acousti c properties in the 



Khomenko et al., 2009; Parchevsky & Kosovichev, 200' 



sphere and interior ([ Hasan e t all 120051: LH anasog e et al 



20071 lHasan fc van Ballegooiienl . 120081 : iFedun et all 12009: 



Shelvag et all 120091 : ISteinerl . 120091 : IVigeesh et all 1200* 
However, due to the non- locality of radiative processes in 
the solar atmosphere, a direct comparison of the plasma 
parameters at a certain geometrical depth in the com- 
putational box with the solar radiation parameters at a 
given optical depth may not be entirely correct. The non- 
locality of radiative transport must be taken into account. 
Khomenko fc Colladosl (|2009h have recently suggested that 
the changes in the height of continuum formation with re- 
spect to the equipartition layer (the layer where the Alfven 
speed is equal to the sound speed, va = c s ), may help to ex- 
plain the appearance of the high-frequency acoustic haloes 
around sunspots. 



Numerical simulations of solar wave phenomena require 
a static magnetic configuration model which incorporates as 
many physical properties of the real Sun as possible. In this 
paper we provide a recipe to create such a model, based on 
the results of numerical modelling of magneto-convection 
in the photosphere. We demonstrate that the spectropolari- 
metric properties of the magnetic configuration we created 
successfully reproduces those of MBPs. In Section 2 we de- 
scribe the technique used to reconstruct the magnetic and 
thermal parameters of the MBP model. The results of the 
spectropolarimetric simulations using the model are pre- 
sented in Section 3. In Section 4 we show our preliminary 
results on the wave mode conversion in the MBP and dis- 
cuss the possible observational signatures. Our concluding 
remarks are presented in Section 5. 
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Fig. 1. Vertical magnetic field component B§ z as a function 
of depth for the photospheric MBPs. A depth level of z = 
corresponds to log T4300 = in the quiet solar model. 



2. Static bright point model 

We use a snapshot from the "plage" magneto-convection 
simulation of t he solar photosphe re undertaken with the 
MURaM code (|Vogler et al.L l2QQ5h to produce the average 
MBP model. Since the average magnetic field flux in this 
snapshot is relatively high (200 G), a large number of inter- 
granular magnetic field concentrations are generated. The 
magnetic field concentrations appear bright in the contin- 
uum (A = 4300 A), and in the G band. This snapshot has 
been used to demonst rate the magnetic n a ture of the G- 
band bright points (Schussl er et al.L 120031 : IShelvag et al.L 
12004 ). To reveal the basic structure of the vertical compo- 
nent of the magnetic field Bq z in the MBPs, we average the 
depth dependences of the vertical magnetic field strength 
over the bright points, which are selected by their enhanced 
G-band intensity and magnetic field strength. Fig. [1] shows 
Bo z as a function of depth. Note that z = on the depth 
scale corresponds to the average optical depth log T4300 = 
for t he quiet Sun, and thus does not include Wilson depres- 
sion ([Wilson fc MaskelvneL Il774f ) . The Wilson depression is 
about 100 km, as will be shown later. 

We use the self-s imil arity assumption 
(Schluter & Temesvarv, 1958; Schussler & Rempel, 



120051 : iGordovskvv fc JainL 120071 :' IShelvag et all l2009h to 
reconstruct the horizontal component of the magnetic 
field. Below we briefly discuss its governing equations. 

For the prescribed vertical magnetic field component 
Bq z (2), defined along the axis of the symmetric magnetic 
field configuration, the (divergence-free) two-dimensional 
magnetic field can be defined as: 



B z (x,z) = B 0z f(xB 0z ) 



and 



B x (x, z) 



dB 0z 
dz 



xf (xB 0z ) , 



(1) 



(2) 



where / is an arbitrary function which describes how the 
vertical component of the field expands with height. 

For the simulations, we choose the following grid pa- 
rameters of the numerical domain in which the magnetic 
field is embedded: the vertical and horizontal extents of 
the domain are 1.4 Mm and 2 Mm, respectively, which are 
resolved in 100 x 200 grid cells. 



In Fig. [2] we show the magnetic field configuration cal- 
culated for this domain using the above definition. The 
function / takes zero values on the side boundaries of the 
domain and is chosen in such a way that the radius of the 
magnetic field region is approximately 200 km at the height 
corresponding to the quiet Sun continuum formation level. 

The magnetic field changes the thermodynamic param- 
eters of the plasma where it is embedded. These changes 
can be quite significan t for the plasma dyn amics and wave 
propagation (see e.g. IShelvag et al.l (|2009h ). To construct 
the magneto-hydrostatic configuration, we substitute the 
values of the magnetic field components obtained using 
Eqs. (pQ) and (j2j) into the equation of magneto-hydrostatic 
equilibrium 



(B- V)B + V 



B 2 



Vp pg. 



(3) 



Here the magnetic field strength is normalised by the factor 
\^4tt. If the magnetic field is prescribed, Eq. (|3]) splits into a 
system of first-order partial differential equations describ- 
ing the horizontal and vertical pressure and density vari- 
ations from their equilibrium values. These equations are 
then integrated numerically. This procedure of integration 
gives satisfactory results in terms of numerical precision. 

We use th e standard model of the solar photosphere 
([SpruitL Il974l ) as the unperturbed background model. The 
unperturbed pressure is recalculated from the standard 
density profile using the hydrostatic equilibrium condition 
Vp M = p u g. Then the perturbations to the pressure p and 
density p, obtained from the solution of Eq. (j3]), are added 
to the unperturbed density and pressure dependencies. 

Partial ionisation effects are taken into account when 
the plasma internal energy and temperature are calculated. 
The system of Saha equations is solved for the eleven most 
abundant elements in the solar photosphere to produce the 
tabulated functions of internal energy e = e(p,p) and tem- 
perature T = T(p,p). Values of internal energy and tem- 
perature are then obtained by interpolation. The density, 
pressure, internal energy and temperature structures are 
shown in Fig. [3j where the lower right panel shows a sig- 
nificant temperature increase above the MBP. This tem- 
perature increase is caused by the increased magnetic ten- 
sion in the upper layers. A similar, but somewhat weaker, 
temperature increase is observed in the dynamic radiative 
magneto-convection simulations. 

For a direct comparison of the radiative properties of 
the model with observations, we need to know the position 
of the continuum formation at some wavelength relative to 
the equipartition layer (va = c s ). The thermal and mag- 
netic structure of the average MBP is shown in the lower 
right panel of Fig. [3j along with the equipartition layer and 
continuum formation level logT43oo = 0. The latter curve 
demonstrates the presence of Wilson depression of the or- 
der of 100 km in the magnetised region, similar to that 
observed for sunspots. The Wilson depression in the MBP 
model is significantly less than that of sunspots (Solanki, 
1993; Watson et al.L l2009f ). The equipartition layer, as is 
evident from the figure, is located deeper than the contin- 
uum formation layer in the region of the strongest magnetic 
field. This means that in the MBP centre, the 4300A contin- 
uum is originating from strongly magnetised plasma with 
va > c s and /3 < 1, and a variety of MHD effects, including 
the magneto-acoustic mode conversion, may be observable 
at this wavelength. 
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Fig. 2. Vertical (left) and horizontal (right) components of the reconstructed magnetic field for the MBP model. 
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Fig. 3. The theoretical density (top left), pressure (top right), internal energy per unit volume (bottom left) and tem- 
perature (bottom right) structures of the MBP. In the bottom right panel, the white line shows the 4300A continuum 
formation level logT43oo = 0, the green line represents the equipartition layer va = c s , and the blue lines are the magnetic 
field lines. 



3. Radiative diagnostics 

The line profile s imulation code STOPRO ([Solankil [l98 



I: 



Frutkerl. l2QQd : iBerdvudna et all l2QQ3t IShelvag et al 

I04, 2007) is used to perform the line profile calculations. 
This code is designed to compute wavelength-dependent in- 
tensities and normalised full Stokes vectors for the atomic 
and molecular line profiles in the LTE approximation. We 
use STOPRO to calculate the continuum intensities at 
4300A and 5000A, as well as G-band intensities for each of 
the 200 vertical rays of the MBP model described in Section 
2. To calculate these we employ a p r ocedu re similar to that 
described in detail in lShelvag et al.l (|2004f ). The G-band in- 
tensities are derived by convolving the simulated spectrum 
with a 10A filter, centered at 4305A, and integrating this 
over the range 4295A - 4315A. The intensities are shown in 
Fig. [U and have been normalised by the values taken in the 
non-magnetic part of the domain corresponding to granula- 
tion. As is evident from the figure, the MBP model appears 
to be bright in its middle x = 0. Darkenings, corresponding 
to the integranular lanes surrounding the MBPs, are found 
at x = —0.3 Mm and x = 0.3 Mm. 

The 6302. 49 A Fe I transition is often used as a diag- 
nostic for the photospheric magnetic field, and the Stokes-I 
and -V profiles for this line are shown in Fig. [5] for two posi- 
tions in the model. The I profile is split by strong magnetic 
field in the centre of the MBP, but the depth of the profile 
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Fig. 4. G-band intensities for the average MBP model. The 
darkenings, corresponding to the integranular lanes, are 
found at x = —0.3 Mm and x = 0.3 Mm. 



is reduced due to a smaller temperature gradient. In the 
weaker magnetic field (Fig. [5j right column), the absorp- 
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Fig. 5. Stokes-I (top) and -V (bottom) profiles calculated 
for magnetic flux concentration at x = 0, (left), and x = 
250 km, (right). 



tion is stronger, and the Stokes- V profile has a larger am- 
plitude. However, the wavelength separation between the 
left and right lobes of Stokes-V remain approximately the 
same in both cases. It should be pointed out that the ob- 
servational profiles are also influenced by the Doppler shift, 
caused by the bulk motions of the convecting plasma. Our 
model does not account for convective velocities and thus 
cannot reproduce the exact central wavelength and C-shape 
of the observed absorption line profiles. 



4. Wave propagation in the MBP 

We use the code SAC to perform simulation s of wave prop- 
agati o n and mode conver sion in the MBP (jShelvag et all 
200$ iFedun et all l2QQ9h . The code solves the full ideal 
MHD equations in a two- or three-dimensional Cartesian 
grid, and has previously been used for a variety of problems 
in the physics of solar oscillations. The boundary conditions 
implemented in the code are "transparent" conditions of the 
simplest type. These conditions lead to an artificial reflec- 
tion from the boundary layers, however, for well-resolved 
perturbations the amplitudes of reflected waves are quite 
small. The full description of the methods u sed, boundary 
cond itions and tests are presented in detail bv lShelvag et al.l 
(2008). For the simulations, the code has been modified 
to incorporate the equation of state and account for par- 
tial ionisation processes. At each time step, the local in- 
ternal energy- and density-dependent adiabatic index I\ 
is calculated using the pre-calculated tabulated function 
Ti = Ti (p, e). The adiabatic index is then used in the equa- 
tion of state to relate the kinetic pressure to the internal 
energy values, thus closing the system of MHD equations. 

Oscillations are generated using a single temporally con- 
tinuous acoustic source located in the non-magnetic part of 
the numerical domain. This source is coherent, acts in the 
horizontal direction v x , extents vertically over the whole 
height of the domain, and its period is T = 30 s. The am- 
plitude of the source is 100 m s _1 , is constant over the 
domain depth and is chosen to be sufficiently small to keep 
the oscillations in the linear regime. The diffusive damping 
introduced into the code by increasing the hyperdiffusion 
coefficients to 0.2, combined with the low source amplitude, 
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Fig. 6. Vertical velocity component in the domain. Lines 
as in lower right panel of Fig. [3l A complicated velocity 
pattern, produced by the interaction of the oncoming wave 
with the magnetic field concentration, is visible above the 
equipartition layer va = c 3 . 



prevent the convectively unstable domain from going into 
the convective regime for a sufficiently long time. 

A snapshot of the vertical velocity component, taken 
at t = 160 s from the start of the simulation, is shown in 
Fig. [6j An oscillatory pattern, bounded by the equiparti- 
tion level va = c s (marked by the green line in the figure), 
is easily distinguishable. This pattern is produced by an 
interference of the slow and fast magneto-acoustic waves, 
propagating differently in the region. Most of the velocity 
pattern is above the continuum formation level logr = 
(white dashed line), and should therefore affect the absorp- 
tion line profiles. 

Fig. [71 shows the 5000A continuum intensity variations 
in non-magnetic media outside the MBP (dashed lines) and 
in the centre of the MBP (solid lines). It is evident from 
the figure that both the absolute intensity variation (upper 
plot) and the intensity variation normalised to the average 
continuum intensity (lower plot), are larger in the MBP 
than in non-magnetic media. This is caused by the larger 
continuum intensity and stronger variation of temperature 
in the deeper (due to depression) layers of the solar pho- 
tosphere in the MBP. Although the source amplitude is 
sufficiently low to keep the linear character of the intensity 
oscillations in the non-magnetic part of the domain (dashed 
lines), a non-linear character of the intensity oscillations in 
the MBP is also observed (solid lines). The steepening of 
the wave front observed in the figure is caused by conver- 
sion of the linear wave into the non-linear regime in the 
partially evacuated plasma of the MBP. 

An oscillatory behaviour can also be detected in the 
magnetic field. To reveal the observational consequences of 
the oscillations in the magnetic field, we apply the standard 
observational technique of measuring the Stokes-V profile. 
We choose a narrow bandpass filter with a width of 0.02 A, 
centered — 0.0525A from the Fe I line centre. We then apply 
this filter to the Stokes-V profiles data, calculated for the 
whole x — t domain. The resultant Stokes-V intensities are 
shown in Fig. [U The vertical component of the magnetic 
field in the region of the line formation is proportional to 
the intensity of the Stokes-V filter. However, it is worth 
noting that the Stokes-V profile for the 6302A line satu- 
rates in regions of strong magnetic field. As Fig. [8] demon- 
strates, the intensity at x = is equal to the intensity at 
about ±0.3 Mm. Also, the magnetic field strength values, 
restored from the Stokes-V intensities at these positions, 
are virtually indistinguishable, thus making the 6302. 49 A 
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Fig. 7. Absolute (upper plot) and relative (lower plot) con- 
tinuum intensity oscillations in the non-magnetic media 
(x = 0.66 Mm, dashed line) and in the bright point centre 
(x = 0, solid line). 



Fig. 9. Stokes-V asymmetries, calculated for the simulated 
time series. 

produced by the perturbation is large enough to be detected 
by modern instruments. 




Fig. 8. Filtered Stokes-V intensities, calculated for the sim- 
ulated time series. 



Fe I line unsuitable for measurements of strong magnetic 
fields. The variation in the filtered Stokes V amplitudes in 
the centre of MBP is about 25-30% in our simulations. This 
value is quite large and should therefore be detectable by 
observations. 

The magnetic field configuration in the simulations 
is localised and non-uniform, and the line-of-sight there- 
fore crosses regions of both zero and non-zero mag- 
netic field. The initial static model does not produce 
an asymmetry in the Stokes-V profiles, since it has 
no velocities w hich are necessary to produce an asym - 
metry (see e.g. | Sankarasubramanian fc Rimmeld (|2002l ); 
iBellot Rubio et all ([2000)). In the presence of a flow, caused 
by the perturbation introduced into the numerical box, this 
magnetic f ield configurat ion can lead to asymmetric Stokes- 
V profiles ( Sol ankilll993f ). The simplest measure of the area 
asymmetry is the integral of the Stokes-V profile over the 
wavelength A = JVdX. The V profile area asymmetry map 
is shown in Fig. [9j The normalised amplitude of Stokes-V 
area asymmetry is of the order of 2% in the center of the 
MBP. This value is close to the observe d values of Stokes-V 
area asymmetries, which are about 3% ([Bellot Rubio et all 
2000j). Thus, the variation of the Stokes-V area asymmetry 



5. Concluding remarks 

We have presented the recipe for creating a static numer- 
ical model of a photospheric MBP. This model allows the 
study of oscillatory properties of such small magnetic con- 
figurations in the solar photosphere directly by measuring 
their radiation intensity and polarimetric properties. The 
magnetic field is extracted from dynamic simulations of so- 
lar radiative magneto-convection and reconstructed using 
the assumption of self-similarity. Thus, our static model in- 
herits many observational properties of the full dynamic 
simulation. 

The model reproduces the G-band brightening in the 
MBP centre. Stokes-I and -V profiles for the 6302A Fel line 
are comparable to solar observations in both the magne- 
tised region and in the region corresponding to the non- 
magnetic solar granulation. The continuum formation level 
in the centre of the MBP is located in the region of strongly 
magnetised plasma, where the Alfven speed is greater than 
that of sound. 

We have used the model to examine the observational 
consequences of sound wave propagation though the mag- 
netic field concentration of the MBP. This preliminary in- 
vestigation shows that the variation of continuum intensity 
is more pronounced in the MBP compared to the average 
granule. Using the radiative diagnostics with the Stokes- 
V profile of the 6302A Fe I absorption line, we demon- 
strate the detectability of the magnetic field variation in 
the bright point and show the appearance of the Stokes-V 
profile asymmetry caused by the oscillations. 

The radiative heating term is not included into the sys- 
tem of MHD equations used to perform the wave dynamics 
study. The absence of the heating term may result in re- 
duced vertical radiative flux in the magnetic flux tube. As 
a result, we have somewhat smaller continuum and G-band 
intensity in the MBPs, compared to the full simulations 
and the real Sun. However, for wave dynamics the effect 
of the radiative term i s expected to be small . As i t has 
recently been shown by lYelles Chaouche et al.l (|2QQ9t l. the 
magnetic flux tubes in the solar photosphere are "reason- 
ably well reproduced" by a thin flux tube approximation 
with no radiative term included. 
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Future investigations will examine the response to 
sources at different locations and with a range of frequen- 
cies. The behaviour and observational signatures of Alfven 
waves in small magnetic elements is particularly important 
when used as a diagnostic for solar plasma parameters, and 
for understanding the energy transport to the corona. Thus 
we plan to extend the model to three dimensions. 
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